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, Abstract 

^ ■ We investigate the integrable higher spin XXZ ehain at the Razumov-Stroganov point. 

We present a method to evaluate the exact value of the eigenvalue which is conjectured 
to correspond to the groundstate of the Hamiltonian for finite size chain from the Baxter 
Q operator. This allows us to examine the exact total energy difference between different 
number of total sites, from which wc find strong evidence for the absence of finite size 
correction to the groundstate energy. 
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1 Introduction 



-,— , The spin 1/2 Heisenberg XXZ chain under the periodic boundary condition 

H{ri) = -- 5]kJcTj+i + ajaj+i + chr?(a|a|+i - 1)]. (1) 
i=i 

is the first discovered and is one of the most fundamental models in quantum integrable 
systems. The Hamiltonian can be diagonalized by the Bethe ansatz [1] to give the eigenvalues 

M 

£:(r/) = J^(2chr?-zz;, -t/;-^), (2) 
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under the following constraints between the parameters {wj} 



nj _ WjWk - 2ch{r])wj + 1 1 
wjWk-2ch{r])wk + lj' ^' 



which is called the Bethe ansatz equation. The groundstate energy can be calculated exactly 
in the thermodynamic limit by making Fourier transform, but it is in general impossible to 
find the exact value for finite size chain, due to finite size correction. However, taking the 
anisotropy parameter to a special value 77 = — 27ri/3, the total number of sites to be odd 
M = 2N + 1 and the total spin to be 5*°* = 1/2 {p = N), the groundstate energy is shown 
to have no finite size correction [2-4] 

£{-27ri/3) = -M. (4) 

For the last ten years, many exact results such as the components of the groundstate wave- 
function and the emptiness formation probability were conjectured to be related to combi- 
natorial objects such as the alternating sign matrix [3-6], some of them have been proved 
[7—9]. Nowadays, the special point of the anisotropy parameter rj = — 27ri/3 is called the 
combinatorial point or the Razumov-Stroganov point. 

In this paper, we describe a method to evaluate the eigenvalue which is conjectured to 
correspond to the groundstate of the Hamiltonian of the integrable higher half-integer spin- 
s XXZ chain of finite size at the Razumov-Stroganov point. By appropriately normalizing 
the overall factor of the Hamiltonian which does not depend on the length of the chain, the 
eigenvalue of the Hamiltonian can be calculated as 

M 

£{r,) = J](2ch(2s7y) - wj - wj'), (5) 
i=i 

under the Bethe ansatz equation 



T-T |_ sh{r])wjWk - sh((2s + l)r])wj + sh((2s - l)r])wk + shij 1 
1 sh{r])wjWk — sh((2s + l)r])wk + sh((2s — l)r])wj + shrj J 

k'^j 



As mentioned before, one cannnot evaluate the eigenvalue exactly for finite size chain in 
general. However, for odd number of sites and at the Razumov-Stroganov point, we recently 
calculated a polynomial called the Q operator whose zeros give the Bethe roots [10]. This 
was obtained by solving the Baxter TQ equation [11, 12] which is essentially equivalent to 
the Bethe ansatz equation. The computed Q operator is based on the variables zj = {wj — 
g-2s»?^y'^g-2s»7y^_^. _ 2^ rather than the variables wj. To evaluate the exact eigenvalue of the 
Hamiltonian from the obtained Q operator, we switch from the Q operator in the variables Zj 
to the one in the variables wj. This leads to the transformation of the symmetric polynomials 
of the variables zj to those of the variables wj, and the simplest symmetric polynomial 
Yl^=i essentially leads to the eigenvalue of the Hamiltonian. By this procedure, we can 
easily evaluate the groundstate energy exactly for finite chain, and one finds strong evidence 
for the absence of finite size correction for higher half- integer spin XXZ chain at the Razumov- 
Stroganov point, generalizing the behavior observed in the spin 1/2 XXZ chain. Calculating 
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the groundstate eigenvalue naively by solving the zeros of the Q operator just only gives 
a numerical value. The transformation of the Q operators is essential to extract the exact 
value, from which one can check the validity of the vanishing of finite size correction. Based 
on the conjecture for the absence of finite size correction, it is enough to consider only M = 3 
sites {N = 1) and M = 5 sites {N = 2) to extract the groundstate energy which can be easily 
evaluated by the Q operator. 



2 TQ equation and Q operator 

In this section, we give a brief review on the results for the Q operator of higher spin 
XXZ chain at the Razumov-Stroganov point [10]. In terms of the variables Zj = (wj — 
e~'^'^'^)/{e~'^^^Wj — 1), the Bethe ansatz equation (6) can be expressed as 



k=l 



or, in the familiar form 



M p 



ysll{Uj - ST]) J ^J^sh{uj - Uk- v)' 

in terms of the variables uj related to Zj by Zj = exp{2uj). 

The Bethe ansatz equation can be obtained from the Baxter's TQ equation 

T{u)Q{u) = sh^'\u + srj)Q{u - r?) + sh^^(n - sr])Q{u + ??) = 0, (9) 

where T{u) is the transfer matrix whose auxiliary space has spin 1/2 and the quantum space is 
the M-fold tensor product of spin s spaces, u is the spectral parameter and rj is the anisotropy 
parameter associated with the XXZ chain. The Q operator Q{u) = 11^=1 sh(u — Uj) encodes 
the information of the eigenstate of the model since the TQ equation (9) reduces to the Bethe 
ansatz equation of the higher spin XXZ chain by setting the spectral parameter u to u = Uj . 
Solving the Bethe ansatz equation is equivalent to computing the Q operator. 

For the half-integer spin s = (L — 2)/2 (L = 3, 5, 7, • • • ) XXZ chain with odd number 
of total sites M = 2N + 1 (N = 1,2,3, ••• ), the transfer matrix eigenvalue of T(n) which 
is conjectured to correspond to the groundstate was found to have a simple form at the 
Razumov-Stroganov point rj = —{L — \)iTi/L. In the sector withp = A^+(2A^+1)(L — 3)/2+m 
Bethe roots, an exact transfer matrix eigenvalue has the following simple form [12] 



TM = 2chf'^-^'<' Vl"V (10) 



2L 

To analyze the corresponding Q operator, it is useful to change to the spectral variables 
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exp(2ti) and Zj = exp{2uj) and redefine the Q operator as 

Q{z) = f[{z - z,) = j^i-iyz^-^ef, (11) 

j=l j=0 

- Yl • • • (j = 1, 2, • • • eff) = 1. (12) 

l<il<i2<-<ij<p 
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We recently calculated the Q operator in two ways. Let us briefly present the results below. 
We restrict to the case m = 0, i.e., in the sector 5*°* = 1/2 [p = N + {2N + 1)(L - 3)/2). 

First, we showed that solving the TQ equation reduces to solving the following set of linear 

(p) 

equations of the elementary symmetric polynomials of Zj 

mm(7V+{27V+l){L-3)/2,^) 

y: n^+'yf+i^^+^)i^--m)=o, (13) 

j=max(0/-2Af-l) ^ ^ 

for£ = 0,l,-- - ,NL + {L-l)/2 {i ^ Lk, Lk + {L - l)/2 {k = 0,l,--- ,N)). 

We also evaluated the Q operator in another way by use of the interpolation formula to find 



LN+{L~l)/2~Lk _ ^Lk\ 



k=0 ^ ' j=0 



for A'^ even and 

{N-l)/2 N 



LN+{L-l)/2-Lk _ Lk\ 



k=0 ^ ' j=0 



, ^ n (.L-l)/2 + Lj , LN-Lk LkML-DnA 



for N odd. 



3 Transformation of Q operators 

It was useful to change to the variables z and Zj to calculate the Q operator. On the other 
hand, it is essential to use the variables w = {ze^^"^ — l)/{z — e^^^) and Wj = (zje^"'' — 1)/ {^j — 
e^^^) to evaluate the exact eigenvalue of the Hamiltonian [4]. This is because E^^^ = X]^=i 
which forms a part of the eigenvalue (5) is the coefficient of wP~^ of the following Q operator 
in the w variables 

X{w) = \{{w - w,) = j2{-iyw^~^Ef\ (16) 

j=l j=0 

Ef = • • • U = 1, 2, • • • ,p), EjP^ = 1. (17) 

l<ii<i2<---<«j<p 
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(p) 

To evaluate the exact eigenvalue, we need to express ' in terms of the symmetric polyno- 

the relation 

(e-4Wi - l)(z - zj) 



(p) 

mials of the variables zj . We insert the relation 



(18) 



into (16) and make the following transformation to express the coefficients of the powers of 
iofei^) 

/ e-4"A-i y Q{z) 

yz_Q-2ni/L j Q(g-27ri/L) 



(p) 

w in terms of e\. 



Q(e-27ri/-^) n ^ _ e-2vri/L 



Q(e- 


-2m/L^ 




1 


Q(e- 






1 


Q(e- 


-2m/L^ 




1 


Q(e- 





J 

p 

Ylie-^'^'/^w -1-Zj{w- e-2-^/^)) 
i=i 

fc=0 
p p—k k 



^ fc=Of=0 j=o ^ ^ ^ ^-^^ 

p p mm{fc,p-a) , , . , 

^ (^_l^c+k^-27Ti{k+p-a-2j)/L( P U^K„P-°^(P) 

fc=0 «=0 j=max(0,fc-a) 



yP-a-jJ Vj, 



(19) 

Equating the coefficients of the powers of w of (16) and (19), one gets 

p mm{k,p-a) . , , , 

E^v) = —l-—y y ^_iy^-2mik+p-a^2j)/L P k W/C\ 

K=U j=max(0,fe— a) 

In particular, we have 

p min(fc,p~l) , 

< = 0?FWriE E (-l)'e--<'«-'-W^(/;* )(:)»''-'el". (21) 

We can furthermore simplify the numerator and the denominator utilizing the relation 

^ ^-iYe% (22) 
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This relation can be shown by combining 



(p)| _ -1 -1 -1 



E 

l<ii<i2<--<jfc<p 

E_ (P) 
Zii ■ ■ ■ Zif_ — , 

l<il<i2<-<ik<P 

which foUows from the z -H- z~^, Zj -H- z~^ invariance of the TQ equation 

N+{2N+l)(L-3)/2 



2ch 



i=i 

Ar+(2iV+l)(L-3)/2 
_^g(l-L)7ri/2L _ g27ri/L ^2Ar+l -|-|- ^ ^ _ ^27rVL 

i=i 

Af+(2Ar+l)(L-3)/2 
_^g(L-l)7ri/2L^^_g-27ri/i^|27V+l "Q 



and 



Q(o)=(-irn^i=i' 

which follows by substituting z = in (14) and (15), as 



E 

l<il<i2<---<ik<P 

n X] ^ii "' \-k 

j=l l<h<i2<--<ip-k<P 

l<n<i2<-<«p-fc<p 



Using this relation, the numerator can be expressed as 

p min(fc,p— 1) 



fc=0 j=max(0,fc-l) 
p-1 



k=0 

and the denominator as 



2e-^V^ 5](-l)^(p - A:)cos ) e(^) 



g(e-2-A) = e-f^/^ ( ^^^^^ ) 



A;=0 



The denominator can also be evaluated by setting z = in (14) or (15). The factor 

g-7rpi/L ^27) and (28) cancels out and we have 

2E ^4(-l)np-fc)cos(^^g^ )ei^) 
ELo(-l)^cos(^)4 

This is the expression of the term ^J^^ = Yf,= 

_^ Wj in terms of the symmetric polynomials 
of the variables z^- The other term X]j'=i ""^7^ which consists another part of the eigenvalue 
of the Hamiltonian (5) gives the same value with 

j2^,=iZ<- (30) 

This can be shown by utilizing the invariance (23) and noting that (29) replaced by wj — t- wJ^ 

and Zj — ;> z~^ holds by comparing ■Wj = (z^e^'^'' — 1)7(2;^ — e^'^'') and w"^ = (^""'^e^^'' — l)/(z~"^ — 



k I- ^--1 



2Ellli-mp-k)cos(--^2h^)e(^) 



ELo(-l)^cos(^)er) 

p 

Y^roj. (31) 



4 Absence of finite size correction 

The procedure to calculate the exact eigenvalue of the Hamiltonian for finite size chain can 
be summarized as: 



(i) Solve the TQ equation to calculate the Q operator in the z variables. This reduces to 
solving the linear equations (13) of the symmetric polynomials e^-^^ of the z variables, or to 
solve the TQ equation explicitly by use of the interpolation formula to get (14) or (15). 

(ii) Relate the Q operator in the z variables and the w variables and express the symmet- 

(p) (p) 
ric polynomials Ej of the w variables in terms of the symmetric polynomials ej: of the z 

variables. 

(iii) Evaluate the exact eigenvalue (5) by inserting the exact value of e^^^ computed in (i) 
to the expression (29) relating E^^^ and e^^^ which is done in (ii). 
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For spin 1/2 (L = 3), it is shown [4] that 

^ 1 N 

by identifying the symmetric polynomials in the w variables with the refined enumuration 
of the alternating sign matrix. (32) can also be observed by the procedure described above. 
This leads to the absence of finite size correction to the groundstate energy at the Razumov- 
Stroganov point 

£{-2TTi/3) = -M. (33) 

Next, we examine the spin 3/2 (L = 5) chain. By calculating e[^^ for various number of 
total sites, we conjecture that the following relation 



holds. This leads to the following conjecture for the eigenvalue 

f(-47rV5) = -^^M, (35) 

which means there is no finite size correction to the groundstate energy for the spin 3/2 chain 

at the Razumov-Stroganov point rj = — 47ri/5 as well. 

(p) 

By examining larger spins, we conjecture that Ej^ can be expressed as 

Y^ujj = A+\2A + cosl^]\N, (36) 



where A depends only on the spin value. One can see ^ = 1/2 for spin 1/2 and A = {l+V5)/2 
for spin 3/2 as above. We have checked this for various spins and total number of sites. This 
leads to the absence of finite size correction to the groundstate energy 

£{-{L - l)Tri/L) = (^{L - 3)cos - 2A^ M. (37) 

Based on the conjecture for the absence of finite size correction, E^^ = X^^^i Wj and the 
groundstate energy can be evaluated by only treating M = 3 sites (N = 1) and M = 5 sites 
{N = 2). Here is a list of Ej=i ^^i ^^i spin 5/2 (L = 7), spin 7/2 (L = 9) and spin 9/2 
(L = 11). Note that the value A in (37) is obtained by substituting = in the expression 
for Yl^=i listed below. 

spin 5/2 (L = 7) 

_ 6(-499 + 525cos(7r/7) + 694sin(7r/14) - 900sin(37r/14)) 

f^^'^^ ~ ~ -235 + 290cos(7r/7) + 350sin(7r/14) - 434sin(37r/14) 

,(2_N) -6 + 120cos(7r/7) + 81sin(^/14) - 38sin(37r/14) 
^ ' -2 + 15cos(7r/7) + 12sin(7r/14) - 6sin(37r/14) ' ^ ' 
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spin 7/2 (L = 9) 



p 



7695 + 43820cos(7r/9) - 26108cos(27r/9) - 32210sin(7r/18) 
351 + 2450cos(7r/9) - 1640cos(27r/9) - 1910sin(7r/18) 



+ {2-N) 



-459 + 250cos(7r/9) - 1360cos(27r/9) - 976sm(7r/18) 
-36 + 40cos(7r/9) - 124cos(27r/9) - 100sm(7r/18) 



(39) 



spin 9/2 (L = 11) 



Y^u,, = {N-l)- + {2-N) 



-75 + 860cos(7r/ll) - 207cos(27r/ll) + 575sin(7r/22) - 378sin(37r/22) + 765sin(57r/22) 
^ -9 + 60cos(7r/ll) - 21cos(27r/ll) + 45sin(7r/22) - 30sin(37r/22) + 55sin(57r/22) ' 
B = -8675 + 9780cos(7r/ll) - 16727cos(27r/ll) + 12895sin(7r/22) 

- 15050sin(37r/22) + 10925sin(57r/22) 

C = -376 + 480cos(7r/ll) - 730cos(27r/ll) + 590sin(7r/22) 

- 670sin(37r/22) + 520sin(57r/22). (40) 

5 Conclusion 

In this paper, we described a method to evaluate the groundstate eigenvalue of the integrable 
higher spin XXZ chain at the Razumov-Stroganov point by the Q operator. Since the spectral 
variables of the Bethe roots convenient to evaluate the Q operator and the Hamiltonian are 
different, we made transformation of the Q operators of the two types of spectral variables. 
This enabled us to get the exact groundstate eigenvalue, which gives strong evidence for the 
absence of finite size correction for higher half-integer integrable spin chain, generalizing the 
behavior observed in the spin 1/2 chain. Calculating the groundstate eigenvalue naively by 
solving the zeros of the Q operator just only gives a numerical value, and the transformation 
of the Q operators is essential to extract the exact value. For the spin 1/2 chain, this 
lead to the identification of the symmetric polynomials of the Bethe roots to the refined 
enumuration of alternating sign matrix [4]. It may be possible to prove for higher spins 
by considering alternating sign matrix for higher spins. Alternating sign matrix for spin 
1 has been considered in [13]. Another promising approach is the investigation from the 
supersymmetric point of view [14, 15]. The eigenstate absent from finite size correction should 
correspond to the zero energy groundstate of a Hamiltonian which have super symmetry. For 
higher spins, progresses have been made in [14, 16]. 
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